Symmetry and environment effects on rectification mechanisms in quantum pumps 
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We consider a paradigmatic model of quantum pumps and discuss its rectification properties in 
the framework of a symmetry analysis proposed for ratchet systems. We discuss the role of the 
environment in breaking time-reversal symmetry and the possibility of a finite directed current in 
the Hamiltonian limit of annular systems. 
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Introduction. Recently there has been a good amount 
of experimental and theoretical activity devoted to study 
quantum pumps [H 0- [3- CO- and quantum ratchets 
|f|. The basic underlying idea is the generation of a net 
current as a response to a time-dependent external field 
without a net static bias. The potential applications of 
this effect captures increasing interest within the com- 
munities of condensed matter physics and chemistry. 

The paradigm of a ratchet system is a device with 
broken spacial symmetry affected by a zero-mean time- 
dependent force. An additional ingredient is the cou- 
pling to an environment, which is usually represented by 
reservoirs or some external noise. An important point in 
the investigation of the ratchet effect has been the un- 
derstanding of the role played by the symmetries in the 
rectification properties of the related devices, fn partic- 
ular, a very simple criterion has been proposed in order 
to decide whether a system driven by a time-dependent 
field is able to support a dc-current |7| ■ It could be stated 
as follows: current rectification is not possible in systems 
where the symmetry operations leading to a change in 
the sign of the relevant current, J{t) — > —J(t), leave the 
ensuing equations for its time evolution invariant. This 
criterion is completely general and applies to both adi- 
abatic and non-adiabatic regimes. Its validity has been 
mainly explored in the framework of classical systems. 

Recent advances of material science, have enabled the 
experimental realization of the ratchet effect in quantum 
pumps P, ■ Charge and spin currents have been gen- 
erated as a response to two harmonic potentials with a 
phase lag, which induce out-of-phase oscillations at the 
walls of a quantum dot. Experiments have been per- 
formed in linear arrays where the quantum dot is in con- 
tact to leads. Under these operational conditions these 
devices are actually open quantum systems and time- 
inversion symmetry breaking is introduced in the prob- 
lem not only because of the phase lag between the poten- 
tials but also through the coupling to the environment. 
An interesting alternative setup is obtained by bending 
the structure to form a ring to generate a dc-current along 
its circumference. An important example of this class is a 
ring threaded by a time-dependent magnetic flux. In the 
case of a flux with a linear dependence on time, Bloch os- 
cillations take place due to the induced constant electric 
field and the coupling to reservoirs is essential to rectify 
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FIG. 1: (Color online) Transmission function for a ring with 
N = 20 sites, h — 9,h = 12 barriers of height Et = 1, 
V\ = Vz = 0.5, Wo = 1 and 8 = 7r/2. Solid red (dark gray) 
and dashed black lines correspond respectively to w' — 0.1, 0. 
Thick and thin lines correspond to ilo = 0.05, 0.3, respec- 
tively. An scheme of the setup is also included. 

this current [1,11. 

In this work we consider the setup of Fig. ^ which 
corresponds to a double-barrier structure embedded into 
a ring. The ring is connected to two reservoirs with the 
same chemical potential fj, and two harmonic potentials 
with a phase-lag are applied at the barriers. Our aim is 
to perform a careful analysis of the relevant symmetries 
of the system on the basis of the scheme of Refs. 
proposed for ratchet systems. We investigate the possi- 
bility of a directed current in the limit where the coupling 
to the reservoirs tends to zero as well as the role of the 
environment in the rectification properties of quantum 
pumps. 

Model. We consider a tight-binding chain of N sites 
with constant hopping w and two barriers of height Eh- 
The ends of the chain are connected to reservoirs. The 
chain is closed with a hopping w' along the bond (IN). 
The Hamiltonian for the full system reads 

H = H 1 +H 2 + H c {t)-w Q Y J {A^ x +H.c.) 

fci 

-wo X]( c iV c te + H.c), (1) 
where H± and Hi are free-electron Hamiltonians with 
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degrees of freedom labeled by k\,k 2 , respectively, which 
represent the reservoirs. The latter are coupled to the 
ring through a hopping wq . The Hamiltonian of the ring 
containing the two oscillating barriers reads 
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H c (t) = -w (4 c '+i + H - c -) ~ w'(c[c N + H.c.) + 
1=1 

[E b + V l cos(n t + 6)}clc h + [E b + V 2 cos(n t)]cl<®) 



This model has two interesting limits: for w' = 0, wp 
it corresponds to the linear array studied in Ref. |10| . 
while for wq = 0,w' ^ it corresponds to the ring 
isolated from the reservoirs. In the latter case, the 
spacial coordinates satisfy periodic boundary conditions 
I + N = I. We assume that reservoirs and barriers are 
symmetrically placed defining a mirror line along a di- 
ameter (see scheme of Fig ^| . Their positions satisfy 
l\ = —1% + 1. The relevant symmetry operations to ana- 
lyze are: 



Si 
S 2 



-l + l, 
-t, 



t-^t- s/n 



(3) 



which cause a spacial inversion combined with a shift in 
the time coordinate, and a time inversion, respectively. 

Isolated annular pump. The ring isolated from the 
reservoirs (wq — 0) defines a Hamiltonian, or closed, sys- 
tem. We now show that if He (t) is invariant under <Si or 
1S2, the directed current along the ring vanishes. 

The evolution of a single particle wave function *f?(t) — 
tpi(t) is determined by (we work in units where h = 1): 
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i olM t ) - E £ lM t )^m(t) = o, 



(4) 
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beinge /im (£) = -Wi5 m , ! i + i-wi-i5 mt i-i + 5i tm [5i t i 1 vi(t) + 
Si,i 2 v 2 (t)], being wi = w, I = 1, . . . , N — 1, u>n — w' and 
vi(t) = E b + Vi cos(n t + S), v 2 (t) = E b + V 2 cos(Ooi). 
The ensuing time-dependent current is 



jr l (t)=ewiJm[^t(t)i'i+i(t)}- 



(5) 



Since the applied fields are harmonic with frequency f^o, 
it is verified Jf sol (t + r ) = Jt sol (t), being r = 27r/fi . 
The dc-component of this current is independent of I, due 
to the continuity condition. It is defined as 



jisol 



1 



1 



/V 



jr\t). (6) 



For Vi = V 2 and <5 = 0, 7r, the matrix elements 
of the Hamiltonian £; )m (t) are invariant under S\ and 
S 2 . Hence, applying Si to the equations (@J, it is 
found ipi(t) = ip-i+i(t ~ 8/£lo), while applying 52, it 
is found ipi(t) — %jji(—t). In the first case, we obtain 
(t) — > —Jf° l (t — 5/D.o), while in the second one, 
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jisol ^ _^ _ jisoi^_^_ j n ^- ne ^- w0 caseSj the final conse- 
quence is J = 0. 

In summary, it becomes clear that, in order to have 
jisol q we neec [ geometrical arrangements with bro- 
ken S\ and 1S2 symmetries. In the case of a symmetric 
static setup like the one we are considering, 1S1 can be dy- 
namically broken by applying time-dependent potentials 
with (i) Vi ^ V 2 and/or (ii) S ^ 0,n. Instead, condition 
(ii) must be fulfilled in order to break S 2 . Altogether, 
we conclude that in the Hamiltonian limit, S 7^ 0, n is a 
necessary condition that the pumped two barrier system 
must fulfill in order to support a finite net current. 

The presence of the environment. We now refer to 
the setup of Fig. 2] with wq ^ 0. As mentioned be- 
fore, the linear arrangement is contained in the limiting 
case w' = 0. In the latter limit the dc-current vanishes 
as Wq — > [lfj. Instead, the above symmetry analysis 
suggests that this may be not the case in the ring geom- 
etry (w' 7^ 0) when S 7^ 0, ir. A convenient theoretical 
framework to study transport phenomena in driven open 
quantum systems is provided by non-equilibrium Green 
function formalism (tj, IToj . For reservoirs with the same 
chemical potential the dc-component of the current 
flowing along the ring reads: 



Ji 



dwf(u)Ti(uj), 



(7) 



where we assume zero temperature, hence, f(u>) — 0(/i— 
lu). The transmission function is [lfj: 

T,(w) = — Wl \w \ 2 [ °^poMlm[Gf 1 (t,w)G^ +1 (w,t) 
+G^ N (t^)Gi, +1 (uj,t)}, (8) 

where Gf" m {t, w) is the Fourier transform of the retarded 
Green function, with respect to the difference of time 
t — t', at the time of observation t. The advanced Green 
function is Gf m (uj,t) = [G^ ; (i,w)]* and po(uj) is the 
density of states of the reservoirs. In this geometrical 
arrangement, there is in general a net charge flow between 
ring and reservoirs. Hence, Ti(u) = T(ui), Ji = J for 
I = 1, . . . , N - 1 and T N {uo) = T'(u), J N = J', being 
J' ^ J, which of course verify Kirchoff rules. The exact 
retarded Green function is the solution of the following 
linear set [TfT |: 
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where Si = S,S 2 = 0, while G^ n (ui) is the retarded Green 
function of the equilibrium ring with barriers connected 
to reservoirs without time-dependent voltages. 
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FIG. 2: (Color online) dc- current J as a function of the 
phase lag for V% — 1 .5 V2 and Wq = 0.5. Thin black and 
thick red (dark gray) lines correspond to Vi = 0.2, 0.5, re- 
spectively. Solid, dashed and dot-dashed lines correspond to 
f2o = 0.05,0.1,0.6. Circles correspond to fits with the func- 
tion J = Ao + Ai sin(5). The chemical potential is fi — 0.1 
and to' = 1. Other parameters are as in Fig. 



We now show results for different pumping conditions 
and strengths of coupling to the reservoirs. We consider 
two infinite tight-binding chains of bandwidth W for the 
reservoirs, described by the density of states po(u>) — 
4-^/1 — oj 2 /W 2 Q(W — uS). All energies are expressed in 
units of w. 

First, we would like to point out that the coupling to 
the reservoirs introduce inelastic scattering events and 
the propagation of the wave-packet along the ring looses 
its coherence. Hence, for strong coupling to the reser- 
voirs, the system evolves smoothly from the annular to 
the linear geometry as w' — > and we expect that the 
transport properties of the ring does not significantly dif- 
fer from those of the linear array. This feature is illus- 
trated in Fig. ^ where we show that for Wq = 1 and 
small w' we can almost exactly reproduce the transmis- 
sion function T(u>) of the linear array. 

In what follows, we consider w' = w and a symmet- 
ric static arrangement with barrier positions l\ and l 2 
equidistant to the reservoirs (l\ = —l 2 + 1). We turn 
to show that a finite net current J may be obtained 
for 6 — 0, 7r if the pumping amplitudes are different 
(Vi 7^ V 2 ). For small Vi,V2, eq. © can be solved per- 
turbatively and when this solution is replaced in JSJ for 
I = N/2, it is obtained 

T(u) ~ u>K| 2 poH{[(y) 2 - (y ) 2 ]7iH 
xIm[G° N/2tl (u + n )G°* N/2+1A {cj + Oo) + 

G°N/2,l( UJ - fi o)G°jv/2+l,l( w - n )] 

H — y~ sm K s ) 72 {(J) x 

[\G% /2 ^ + n )\ 2 -\G° N/2jN (u;+Q )\ 2 

-\G° N/2A (w - n Q )\ 2 + \G° N/2 ^{u - n )| 2 ]}, (10) 

being 7i( w ) = \G° itl (u)\ 2 + \G^ n (lu)\ 2 and l2 {uj) = 
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FIG. 3: (Color online) dc- current J as a function of the 
coupling to the reservoirs. Upper panel corresponds to Vi = 
I.5V2 and 5 = and lower panel corresponds to Vi — V2, being 
Vi — 0.5 and S = n/2. Circles, squares, diamonds, upper and 
lower triangles correspond to Qo = 0.01,0.05,0.45,0.6,0.75, 
respectively. Other parameters are as in Fig. [5] 



BB[Gf ul {u)[G^ uN {ui)]*]. Therefore, it is found that 

the net current behaves as J ~ i?o[( 2 ) 2 — (i^) 2 ] + 
B1V1V2 sin(<5), being Bq,Bi real coefficients. The exact 
solution of J as a function of S is shown in Fig. For 
the smallest Vi, there is agreement with the functional 
behavior suggested by (|10fl . For higher V\, departures 
from this behavior are observed. In any case, the feature 
we want to emphasize is that for 5 = 0, it, symmetry Si 
is dynamically broken for V\ ^ V 2 , but £2 is still an exact 
symmetry of the Hamiltonian He (t) for the isolated sys- 
tem. A non-vanishing J at these points is a consequence 
of the fact that the latter symmetry is broken due to the 
coupling to the reservoirs. Another issue worth mention- 
ing is that the functional behavior of J we are finding is 
just the one observed in the experimental work 

The role of the environment in breaking time-reversal 
symmetry is highlighted in Fig. [21 For finite Wq, S 2 is 
broken due to the coupling to the reservoirs. The figure 
illustrates the behavior for Si dynamically broken in two 
different ways: The upper panel corresponds to S = and 
Vi 7^ V%, while the lower one, to V\ — V 2 and 6 — n/2. In 
the first case symmetry 1S2 is restored as uio ^ 0. Thus, 
J — > as the system evolves towards the Hamiltonian 
limit. Instead in the second case, £2 remains broken at 
wq — and J may achieve a finite value. 

Finally, let us comment on the behavior of J as a func- 
tion of the pumping frequency Slo- For low V\,V 2 , an 
expansion of (|10|) in powers of flo, leads to J oc fio. For 
completeness, we also show in Fig. SI the exact behavior 
of J for arbitrary fio, Vi and V2. A large | J| is obtained 
when fio is resonant (i.e when it coincides with the en- 
ergy difference between two levels of the isolated system) . 
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FIG. 4: (Color online) dc-current J as a function of the pump- 
ing frequency Qo for S — tv/2. Thin and thick lines correspond 
to couplings to the reservoirs Wq = 0.5, 1, respectively. Blue 
(black) solid and red (dark grey) dashed lines correspond to 
fi = —1.14 and Vi = V2 = 0.2, 0.5, respectively. Black dashed- 
dotted lines correspond to fi = —1.14 and V2 = 1.5Vi with 
Vi = 0.5. Magenta (gray) two dots-dashed line corresponds 
to Vi = V2 = 0.2 and fi — 0.1. Other parameters are as in 
Fig. |21 The inset shows a zoom for small fio- 

The inset shows that J changes linearly in 11q for small 
enough pumping frequencies, which is typical of adiabatic 
driving |3|,ul ■ Such a linear behavior is, however, not ex- 
pected when the coupling to the environment vanishes. 
This is because in that limit G^ n (a->) corresponds to a 
sequence of poles at the energy levels of the free ring and 
it is not possible to perform a power expansion of (|10ll . 
Furthermore, for small pumping amplitudes, the struc- 
ture of © suggests that J is only sizable for resonant 
fio, in agreement with Ref. 

Discussion. Previous discussions in the literature on 
the role of symmetries in quantum pumps suggest dif- 
ferent behavior in adiabatic and non-adiabatic regimes. 
The idea of adiabatic driving is associated to small pump- 
ing amplitudes and low frequency Qq compared to the 
inverse of the typical time for the particle propagation 
through the device. While in practice this definition im- 
plies J oc fio, this concept is usually formulated in terms 
of some approximation for the scattering matrix 0, Q 
and it has been pointed out that such definition strictly 
applies only to isolated systems 0- 

The symmetry analysis carried out in this work leads 
us to conclude that in open quantum systems the fun- 
damental condition to be fulfilled, in order to obtain a 
dc current, is spacial inversion symmetry breaking in 



the Hamiltonian Hc(t). Since we have not introduced 
any assumption regarding the pumping amplitudes and 
frequencies, we conclude that this condition should ap- 
ply to both adiabatic and non-adiabatic regimes. Let 
us support with examples the fact that this conclusion 
applies, in particular, to the adiabatic regime. We can 
mention, at least, three models with time inversion in- 
variance in the Hamiltonian limit but broken spacial in- 
version symmetry where J behaves linearly in the pump- 
ing frequency: (i) The system considered in the present 
work with 5 = 0, ir, V\ ^ V%, finite wq and arbitrary (even 
zero) w' (see Fig. 0J. (ii) Only one pumping potential 
applied away from a symmetric point under spacial in- 
version. In |10| this problem has been solved in a linear 
array and identical procedure and conclusions apply for 
the annular geometry with contacts to reservoirs, (iii) 
A ring threaded by a linear time-dependent flux coupled 
to reservoirs. This basic problem generated interesting 
discussions some time ago on the nature of resistive be- 
havior More recently it has been exactly solved 
For low enough driving, the dc current is linear in the 
induced emf, which is the effective pumping frequency 
of this problem. In all these systems, the key point is 
that time-inversion symmetry breaking is introduced by 
their coupling to the environment. Another important 
conclusion is that time-inversion symmetry breaking in 
the Hamiltonian Hc(t) is, instead necessary to obtain a 
dc-current in the isolated ring. In order to achieve this, 
a minimum of two time-dependent voltages with a phase 
lag non-commensurate with n are needed. 

To finalize, we mention that the setup of Fig ^ with 
vanishing w' and finite wq can be viewed as a schematic 
model to capture the role of symmetries in the transport 
properties of the system studied in Ref. lj. Our re- 
sults are consistent with the behavior J = A + B sin (5) 
observed in that experiment and suggest that the small 
finite J at 6 = 0, tt can be naturally explained as a con- 
sequence of a slight difference in the amplitudes of the 
pumping voltages. 
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